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O 1 ABSTRACT 

r Jj Context. Gas within a galaxy is forced to establish pressure balance against gravitational forces. The shape of an unperturbed gaseous disk can 

, ^ , ■ be used to constrain dark matter models. 

Aims. We derive the 3-D H i volume density distribution for the Milky Way out to a galactocentric radius of 40 kpc and a height of 20 kpc to 

— < , constrain the Galactic mass distribution. 

Methods. We used the Leiden/ Argentine/Bonn all sky 21 -cm line survey. The transformation from brightness temperatures to densities depends 
on the rotation curve. We explored several models, reflecting different dark matter distributions. Each of these models was set up to solve the 
combined Poisson-Boltzmann equation in a self-consistent way and optimized to reproduce the observed flaring. 

Results. Besides a massive extended halo of M ~ 1.8 10 12 Mq, we find a self-gravitating dark matter disk with M = 2 to 3 10 1 1 Mq, including 
a dark matter ring at 13 < R < 18.5 kpc with M = 2.2 to 2.8 10 10 Mq. The existence of the ring was previously postulated from EGRET data 
■ and coincides with a giant stellar structure that surrounds the Galaxy. The resulting Milky Way rotation curve is flat up to R ~ 27 kpc and 
Cj^> slowly decreases outwards. The Hi gas layer is strongly flaring. The HWHM scale height is 60 pc at R = 4 kpc and increases to ~ 2700 pc at 
r/"*- , R = 40 kpc. Spiral arms cause a noticeable imprint on the gravitational field, at least out to R = 30 kpc. 

' Conclusions. Our mass model supports previous proposals that the giant stellar ring structure is due to a merging dwarf galaxy. The fact that 
the majority of the dark matter in the Milky Way for R £ 40 kpc can be successfully modeled by a self-gravitating isothermal disk raises the 
• <"H , question of whether this massive disk may have been caused by similar merger events in the past. The substructure in the Galactic dark matter 
' disk suggests a dissipative nature for the dark matter disk. 

Key words. Galaxy: disk - Galaxy: structure - Galaxy: kinematics and dynamics - galaxies: interactions - ISM: structure 



1. Introduction 

The shape of the Galactic H i disk is interesting for its own 
sake. Size matters and we would like to know how the Milky 
Way compares to other galaxies. When trying to convert the 
observed brightness temperature distribution to a volume den- 
sity distribution one realizes, however, that this transformation 
depends on the shape of the rotation curve. Determining the ro- 
tation curve needs accurate independent measurements of ve- 
locities and distances for a broad distance range, a daunting 
task. The first large-scale m ap of the H i volum e density dis- 
tribution in the Milky Way ( Westerhoutl 1 1957b was therefore 
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generated from a mass model ( Schmidt , 19561) . For this model, 
the rotation curve at distances R ^ 10 kpc is falling. 

Observations, however, give opposite results. 

' 19761) were the first to notice a 



Georeelin & Georgelin 



give opposite 
were the first to 



rising rotation curve for a relatively large sample. They used 
H a radial velocities and distances of stars that excite Galactic 
Hn regions and found an increase of v ro t for R 10 kpc. 
Additional systematical variations in v ro t for the northern and 
southern parts of the Galaxy led to the question of whether 
a rotation curve from an axisymmetric model suffers from 
systematic errors. 

The discrepancy between mass models and observations 
has persisted until recently. Theo retical mass mode l s pre 



diet a falling rotation curv e (e.g. iBahcall & Soneir; 



I Caldwell & Ostrikeii 1 198 it Ikohlfs & Kreitschman 
Haud & Einastol Il989t iDehnen & Binnevll 1998b while obser- 



1981 



in, 



1980; 



vations have consistently indicated the opposite (e.g. Blitz, 
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1979 


Schneider & Terzian, 


19831 Fich et al. 


1989 


; Merrifield, 


1992 


Brand & Blitz, 1993; 


Honma & Sofue 


1997 


). 



The situation is further complicated by the fact that the ro- 
tation curve depends on the choice of Galactic constants. A 
high rotational velocity at the position of the Sun leads to a 
rising rotation curve, while the opposite holds for a low veloc- 
ity. The interrelation between Galactic constants, dark matter 
models, and the shape of the H i distribution has been studie d 
inaseries of papers bv lblling & Merrifieldldl998l. r2000. 2001). 



These authors constrain the Galactic constants to Rq S5 8 kpc 

an d vrri ^ 200 km s -1 . 

Rohlfs & Kreitschmannl dl988l) considered models with a 



rising rotation curve for 11 ^ R ^ 18 kpc but found it difficult 
to explain this by a dark corona. Such a solution would result 
in a high rotational velocity of 400 S= v ro t ^ 500 km s _1 outside 
the Galactic disk, which appears highly implausible. In addi- 
tion, such a dark matter distribution would disable an exponen- 
tial stenar_djsk 2 _which conflic ts with observations. 



Binney & Dehnenl (119971) discuss the question of whether 
distance errors may mimic a rising rotation curve. They pro- 
pose that most of the tracers that appear to be at R ^ 11 
kpc might actually be concentrated within a ring at R ~ 
14 kpc. The resulting rotation curve would be roughly con - 
stant or gently fall i ng. Several groups -iNewberg et all ( 2002); 
Ibata etal.1 d2003l): lYannvetalJ J2003I); iMartin et al.1 (12005b : 
Martfnez-Del gado et al.l d2005l) : lBellazzini et all d2006l) - have 
subsequently claimed to have detected stellar streams in the 
Milky Way, right at the region of interest. A ring-like mass con- 
centration within this area plays an important role throughout 
our paper. 

The IAU -recommended Galactic constants Rq = 8.5 kpc 
and Vf?) = 220 km s -1 fit best with a flat rotation curve 
dFich et all fl989^ Typical Sb g alaxies have nearly flat outer 
rotation curves (Sofue & Rubin] |200lJ). It seems to be an ac- 
ceptable compromise to assume that the Milky Way is typical 
in this respect. For about a decade now there has been little dis- 
cussion about the Galactic rotation curve. A flat rotation curve 
with v ro t = 220 km s _1 was used for the most recent determina- 



tion of t he Galactic volume density distribution by lLevine et al 
d2006al) . 



The need to obtain better constraints for the radial mass dis- 
tribution in the Milky Way has been stressed by many authors. 
Gas can play a particularly interesting role as a tracer of the 
gravitational potential, provided that its velocity dispersion is 
position-independent and that there are no external forces. The 
inner part of a galactic disk is dominated by stars. This results 
in strong gravitational forces k z perpendicular to the disk. The 
influence of the stars decreases with increasing galactocentric 
distance R. That allows the gas, which withstands the forces 
by means of internal turbulent motions, to reach larger scale 
heights. The H i disks in galaxies are typically three times more 
extended than their stellar counterparts. This property, as well 
as the fact that the outer parts of galaxies are dominated by dark 
matter, highlights the im portance of the gaseous disk as a trace r 
of the mass distribution ( Sofue & RubinT 2001 ; Combesl 12002). 



The Galactic outskirts are very interesting in this respect. 

Gase ous flaring wa s first used as a diagnostic of halo prop- 
erties by Oiling! d 19951) . who showed that the halo of NGC 4244 



was highly flattened. iBecquaert & Com bes d 19971) obtained a 
comparable result for NGC 891. Applying a similar analy- 
sis to the Milky Way, however, has led to conflicting results 
dOlling & Merrifiel(j|l998ll2000ll200lb . lNaravan et all d2005l) 
show that neither an oblate nor a prolate isothermal halo can 
explain the Milky Way flaring. They advocate a steeper fall-off 
for the halo density. 

Surprisingly, little is known about the Galactic gas- and 
mass distribution at large radial distances. Our knowledge is 



limited to distances R ^ 22 kpc dBinney & MerrifielcU Il998 



and references therein). The Galactic outskirts are extremely 
sensitive to the dark matter distribution. The main driver of 
the present paper is therefore to explore the mass distribu- 
tion at large radial distances. A second, even more severe 



problem has been pointed out by Dehnen & Binneyl (|1998). 
There are large observational uncertainties in the mass distri- 
bution perpendicular to the disk. Constraints at z ^ R are re- 
quired. The gravitational acceleration k- in the solar vicinity 
is well known from stellar statistics up to distances of z £ 
1 .5 kpc dKuiiken & Gilmorel 1 199 it iHolmberg & Flvml l2000t 
Korchagin et all 120031: iHolin berg & F lvnnL |2004. but little is 

)f If 



known at larger distances. Figure 1 1 of [Holmberg & Flvnn 
(2004) impressively demonstrates that observational con- 
straints are missing from star counts. Extra-planar Hi gas, 
though fain t, may be a better probe since it extends to larg e 
z distances ( Lockman & Gehmanlll99ltlKalberla et allll998 ). 



This paper is organized in the following way: Sect. 2 de- 
scribes our data analysis, in particular the methods used to de- 
rive the Hi flaring. In Sect. 3 we discuss how mass models 
affect the derived Hi volume density distribution and flaring. 
Each mass model is checked for consistency. The rotation curve 
associated with the best-fit model is discussed in Sect. 4. We 
consider further consistency checks for our model, using flar- 
ing data from the literature and HIPPARCOS surface densities 
in Sects. 5&6. The warp parameters are summarized in Sect. 7. 
We discuss deviations by the Milky Way flaring from the ax- 
isymmetric model and spiral structures in Sect. 8. Uncertainties 
and possible biases due to model assumptions are considered in 
Sect. 9. The results are summarized in Sect. 10. We discuss the 
nature of the dark matter ring and disk in Sect. 1 1 . 



2. Data analysis and methods 



Our study is bas ed on the Lei d en/A rgentine/Bonn (LAB) 
Hi line survey dKalberla et all 120051) . This survey com- 



bines the southern sky su rvey of the Institu to Argentino 



de Radioastronorma (IAR) dBajaja et all 120051) with an im 



proved version of the L eiden/Dwingeloo Survey (LDS) 
dHartmann & Burtonl Il997l) . Currently, this is the most sen- 
sitive Milky Way Hi line survey with the most extensive 
coverage both spatially and kinematic ally. The line profiles 
have been corrected for spurious side-lobe emission. Low- 
level emission features can be interpreted safely as tracers of 
Galactic structure out to large distances, and we intend to ex- 
plore the outskirts of the Milky Way. 
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2.1. T B (l,b,v) ton(R,z,(f>) conversion 

We use cylinder coordinates R, z, <p to describe the Milky Way 
disk; <p = is in direction I = 0. We use the IAU recommen- 
dations for the galactic constants; R0 = 8.5 kpc and V0 = 
220 kms~'. All the individual steps necessary for convert- 
ing the observed brightness temperature distribution Tb{1, b, v) 
as a function of galactic longitude, latitude, and LSR veloc- 
ity to densities n(R, z, 4>) ha v e been described before i n grea t 
detai l, e.g.. bv IWesterhoutj ( 1957 ): iHenderson et ail |l982); 



phase medium. We assumed that the streaming velocity v rot is 
dominated by gravitational forces 



Burton & Te Lintel Hekkert (1986); Diplas & Savage (1991); 
VoskesI dl999l) : INakanishi & Sofuel d2003l) and iLevine et all 
(2006a). We sketch our procedure only briefly; but where we 
deviate from other authors, we discuss those points in detail. 



Previously, ILevine et al.l d2006al) have used the Hanning 
smoothed version of the LAB survey to derive n(R,z, 4>). Our 
approach is similar, but we made use of the original telescope 
data, cleaned for instrument al problems (the vers ion prior to 
the generation of FITS maps; Kal berla et alj d2005l) ). 



To derive volume densities n(R, z, <p) for each volume ele- 
ment centered at Ro,zo,<po we calculated the emission within 
this cell. We first determined which of the observed Tb{1, b, v) 
profiles may show some contribution from this volume ele- 
ment. For each profile we calculated column densities along 
the path length within the cell under consideration. The volume 
density was obtained after division by the path length. The pro - 
cedure is similar to the one used bv iDiplas & Savagd dl99ll) . 
We assumed that the gas is optically thin. This may, however, 
be incorrect for som e parts of the Galactic disk at low latitudes 



(Burton_etaLL|1992). 



The LAB data are expected to predominantly consist of 
emission from the Galactic Hi disk, gas in a roughly cir- 
cular motion. But there are contributions from the local gas, 
from galaxies, also from high- and intermediate velocity clouds 
(HVCs and IVCs, respectively), which are not directly con- 
nected to the disk. Different attempts have been made to ex- 
clude such features from the data base right at the beginning 
of the data processing. Some authors have interpolated and 
subtracted the local H 1 emission. Automated filtering of the 
Tb(1, b, v) data may cause unknown biases. We decided not to 
filter the observational data base but rather the n(R, z, <f>) distri- 
bution. The details for our procedure are given in Sect. 2.2. 

Most important for the conversion of Tb(1, b, v) to n(R, z, <p) 
is the translation of velocities to distances via a rotation curve. 
The exact shape of the Milky Way rotation curve is controver- 
sial, so some assumptions are needed to overcome this deficit. 
Any prejudice at this point may bias the results, so we decided 
to study how different rotation curves, corresponding to differ- 
ent mass models, may affect the results. We also used a gener- 
alized rotation law v mt {R,z), allowing the azimuthal streaming 
of the gas above the disk to deviate from the rotation of the 
disk. 

A strictly circular motion of a particle at a constant z dis- 
tance above the disk is not possible. We considered the proper- 
ties of a gas layer and describe the average rotation for a vol- 
ume element with constant average volume density in a multi- 



, t (R,z) = ^Rd®(R,z)/dR. 



(1) 



Such a relation implies, for constant galactocentric radius R, 
a decreasing rotation velocit y with increasing z distance. It 
was used for the first time by IWesterhoutl d 19571) to derive the 
Galactic H 1 volume density distribution, but later publications 
have assumed cylindrical rotation. 

In comparing cylindrical with lagging rotation we recall the 
cond itions that must apply for the cylindrical case dKalberla , 
2003, Paper I): the density distribution has to be barotropic, 
p = f(p). A barotropic fluid is an idealized fluid in which 
the pressure is a function of only the density. The interstellar 
medium, however, is dominated by several co-existing phases. 
Here we need to consider at least two distinct Hi phases, a 
cold neutral medium, and a warm neutral medium (CNM and 
WNM, respectively). These have ve ry different densitie s but 
are basically in pressure equilibrium dWolfire et all 120031) : and 
phase transitions are possible. Barotropic conditions do not ap- 
ply to such a multi-phase medium, which probably dominates 
most of the volume under consideration. It is therefore highly 
questionable whether the common assumption of cylindrical 
rotation applies for the Hi gas in the Milky Way disk. We 
therefore exercise the option of using both the cylindrical and 
lagging rotations for our reduction. For a more detailed discus- 
sion, see Sect. 4.1. 

Except for replaci ng v m t(R) by Vmt(R, z) , we used the pro- 
cedure as described in ILevine et al.l (|2006a, Eqs. 1&2) includ- 
ing their first order epicyclic streamline correction; which is a 
great improvement over previous distance derivations and im- 
portant for overcoming discontinuities at large radial distances. 
No streamline correction was applied to R < 8.5 kpc, since 
the correction is not applicable there. Inside the solar circle 
the translation between distance and velocity is ambiguous be- 
cause there are two positions along the line of sight that have 
the same velocity. These positions, however, may have different 
z distances, so their densities may differ accordingly. 

We modeled the H 1 distribution, taking the expected flar- 
ing according to the mass model into account. Incrementing 
along the line of sight, we calculated expected densities at near 
and far distances and weighted the derive d quantities accord- 
ingly . Our method is comparable to that of INakanishi & Sofue 
(2003k excep t that they use a sech 1 approach according to 
ISpitzerl (1 19421) in place of a more detailed model. We applied 
no correct ion for densities around t he tangent points, as de- 
scribed by INakanishi & Sofuel (I2003L Sect. 3.2) since this tends 
to cause discontinuities. The focus of our paper is on data at 
R > 8.5 kpc. However, as data inside the solar circle are found 
to be useful for checking the consistency of a model, they are 
included in most of the plots. 

The n(R,z,(p) distribution was calculated and interpreted 
for < R < 40 kpc, -15 < z < 20 kpc at a resolution of 
AR = 100 pc, Az = 100 pc, and A<p = 1°. In addition we al- 
lowed /, b, v to be traced back for each position R, z, <p- To avoid 
biases due to local H 1 features we flagged and discarded all 
data with \b\ > 30°. Similarly we flagged data with |v L sr| < 10 
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kms -1 . In the direction of the Galactic center, we discarded 
data originating from \R\ < 3.5 kpc. This region, approximately 
with the extent of the bar, is affected by ambiguit ies due to large 
velocity dispersions (Wei ner & Sellwoodill999l Fig. 8). 



2.2. Determination of the first moments 

To derive genuine disk gas we initially flagged all data for 
\b\ > 30°. But then we did not discard any features of the 
observed spectra to distinguish between gas within the disk 
and gas that is, most probably, not associated with the disk. 
The derived volume density distribution contains both contri- 
butions. An inspection of the 3-D distribution shows that it is 
relatively easy to distinguish between both contributions visu- 
ally, if one generates video sequences. The disk component, in 
circular motion, stays mostly stationary, while there are other 
components that appear to move toward the Sun as the frames 
approach the Sun or, opposite, as one moves away from the 
Sun. Galaxies, in particular, show up as fast-moving isolated 
spots. This effect is most prominent if the video frames loop in 
azimuth angle <p. 

To focus on the stationary disk emission, we averaged 
n{R,z,4>) over 15° in <p. This average does not have a major 
effect on the stationary disk component but degrades spuri- 
ous emission caused by H i gas in non-circular motion con- 
siderably. We therefore used it as an initial guess of the first 
three moments; the surface density £, the first moment zq, rep- 
resenting mid-plane, and the second moment cr, representing 
the scale height of the gas. As mentioned before, data with 
\b\ > 30°, [vlsrI < 10 km s _1 , had been flagged in the first place. 
In addition, we flagged all data in the original data base outside 
the range zq - 2.35 cr < z < Zo + 2.35 cr. For a Gaussian dis- 
tribution, 2% of the data would be discarded by this criterion. 
On average we find more than double that amount. As a final 
step we determined the moments I.(R, </>), zo(R, </>), and cr(R, <f>) 
from those parts of the unsmoothed data that was not flagged, 
thus excluding extra-planar gas. 

We verified each of the reduction steps by inspecting the 
data cubes visually. Most critical is the flagging process. We 
found the automated procedure worked as expected for a ma- 
jority of cases. However, many features exist that provide no 
clear evidence whether or not they belong to the disk. Figures 
Q] & [2] show example maps for the volume density distribution 
at azimuth <p — 110° (right hand side) and opposite, <p = 290° 
(left hand side). Most of the flagging of extra-planar gas starts 
already in regions below \b\ $ 30°. The latitude limit \b\ = 30°, 
naively introduced to blank out local gas at high latitudes, gets 
justified this way and does not affect our analysis in a quantita- 
tive way. 

Figure Q] displays volume densities for our best-fit model, 
while Fig. [2] was derived from the standard model with a con- 
stant rotation curve; for a more detailed discussion of the dif- 
ferences we refer the reader to Sect. 4. 1 . Prominent in Fig. [T]is 
the elongated feature at 12.5 5 R ^ 16.5 kpc and 4.3 ^ z ^ 6.3 
kpc. This is the HVC complex R located at 62° 55 I 5 73° and 
5°.5 < b < 15° with radial velocities of -156 < v LSR £ -100 
kms 1 . There are some indications that this HVC complex 



may be interacting with the disk ( Kal berla & Haud , 2006), but 
this feature cannot be considered part of the stationary galac- 
tic disk. It was removed by our automatic flagging. Similarly, 
some components showing up in Tb(1, b, v) as intermediate ve- 
locity clouds (IVCs) were removed by the flagging. 

Figures Q] & [2] are examples only, but they do show clearly 
that warp and flaring of the Hi disk are asymmetric. In par- 
ticular, as discussed in Sects. 8 & 9 in more detail, the region 
90° S ^ S 110° appears strongly perturbed for R ^ 15 kpc. 
For this reason we exclude this azimuth range in the following 
whenever we determine average properties of the H i distribu- 
tion. 

Finally, we emphasize that our data reduction is fully auto- 
mated, without any human interaction in the flagging process. 
The procedures were identical for all the mass models and ro- 
tation curves we considered. 

3. Mass models and volume density distributions 

Probably the most frequently used rotation law for the 
Milky Way assumes a con stant rotation v mt (R) = 220 
kms 1 ( Brand & Blitz . 1993b . This naive model is completely 
free of specific assumptions about the dark matter distribution. 
We use it for reference. 

We display in Fig. [2] the volume density distribution de- 
rived from the reference model at azimuth 0=1 10°. The most 
prominent features, considered to belong to the disk, are in 
good agreement with Fig. Q] for our best-fit model. It is the 
changes for HVC complex R that are remarkable. This fea- 
ture has almost disappeared. The main emission from this HVC 
complex has shifted to an azimuth of 104°, while at the same 
time the distance in R has increased by 2 kpc, in z by 1 kpc. 

3.1. Definition of self-consistent mass models 

Our mass model contains the following components: an 
isothermal stellar disk, represented by thin and thick disk com- 
ponents, a bar and a bulge and isothermal gaseous disks for 
molecular gas, the cold neutral medium, the warm neutral 
medium, diffuse ionized gas (DIG), a hot ionized, and a cold 
neutral halo component. Surface densities and scale heights 
for all of these components are adjusted to match the obser- 
vational constraints (we refer to Paper I for details). Such a 
distribution can be supplemented by any dark matter distribu- 
tion that matches known observational constraints. Most im- 
portantly it mus t be consistent w ith the dark matter distribution 
on large scales dZaritskvlll999l) . As a constraint for isothermal 
dark-matter halo models, we demand a total mass of ~ 2 10 12 
Mq within R ^ 350 kpc. We demand V0 = 220 kms'for the 
local dark matter distribution. 

The mass distribution for each model is checked for self- 



consistency as described in Paper I (it has to solve the com- 
bined Poisson-Boltzmann equation). The visible mass distribu- 
tion (stars and gas) is allowed to adjust its scale height to the 
total gravitational potential from disk and halo until conver- 
gence is achieved. The model disk is flat, but for comparison 
with observations, we applied a coordinate transformation to 
take warping into account (Eq. |4j». 
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Fig. 1. Derived volume density distribution for our best-fit model at azimuth 0=1 10° (right hand side) and opposite, <p = 290° 
(left hand side) for < R < 40 kpc and -15 < z < 20 kpc. Densities affected by local emission with IvlsrI < 10 kms -1 are 
flagged black. The curved white, thick lines indicate latitudes \b\ = 30° ± 0.5°. The derived mid plane of the Hi distribution is 
indicated by the thin white line close to z — kpc, and the red dots represent the dispersion of the flaring H i gas layer. The 
outer white dots indicate the initial estimate for the extension of H i gas considered to belong to the disk. A logarithmic transfer 
function was chosen for n < 0.5 ctrT 3 to emphasize low densities. 



SO 



v = 230 kjji/bi, cyliiidrlo ru La I ion 




Fig. 2. Volume density distribution derived under the assumption of a constant cylindrical rotation with v rot = 220. See Fig.Q]for 
further explanation. 



3.2. Self-consistent H i flaring curves 

For each mass model one has to verify whether it is consistent 
with the derived H i distribution. In each case we determined 
the rotation curve, converted Tb(1, b, v) to n(R, z, <p) and calcu- 
lated the moments T.(R, (f>), zo(R, (f>), and cr(R, (/>). We then com- 



pared the average flaring h{R) = f < cr{R, <p) >, derived from 
the data, with the expected flaring. For convenience we used 
for the gaseous flaring / = V2 ln2 to describe the half width 
at half maximum (HWHM) of the H i layer. This convention 
allows an easy comparison with the literature values. 
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3.2.1 . Flaring of the v rot = 220 km s 1 reference model 

Figure [3] shows the flaring curve derived for a flat, v ro t = 220 
kms -1 rotation curve. We excluded 90° ;$</>;$ 110° and those 
regions in direction of the center or anti-center that might be af- 
fected by spurious emission not connected with the H i disk. To 
allow a comparison of our flaring data with previous determi- 
nations we include in Fig.[3]some data from the literature . The 
data for R < R© (O) were derived from lCelnik et al. ( 1979b by 
averaging over ten positions from their Table 3. The literature 
data agree well with our results, except for 7 £ R £ 11 kpc. 
Our analysis is insensitive within this range. 



■■ 220 km/s, cylindrical rotation 



3 - 



2.5 



1.5 - 



0.5 



north 




average 




south 




Wouterloot et al. • 




Merrifield >— *— 




Celnik et al. O 

















10 



15 20 25 
R (kpc) 



30 



35 



40 



Fig. 3. Derived Hi flaring for the standard model. The thick 
red line gives the mean extension of the flaring Hi layer 
(half width at half maximum). The dotted blue and green lines 
give the flaring in the northern and southern hemispheres, re- 
spectively. From t he literature we include : Q derived from 
Celnik et al.ldl979l). X with e rror bars from lMerrifieldl dl992h , 
and« IWouterlootet alJ dl99dL Table 1). 



3.2.2. Flaring caused by a spheroidal halo 

We tested several conventional halo models, some of 



which were already discussed in Paper I. First we dis- 
cuss non-singular spheroids with flattening q (see e.g. 
Oiling & Merrifield 1200 lllNarayan et alll2005h . 



p(R, z) = po 



R 



R 2 C +R 2 + (z/q) 2 



(2) 



Oblate halos have q < 1, prolate q > 1. po is the density at 
the center, R c defines the core radius, and p is the power index 
dNarayan et all 12005 ). We consider the case p — 1 first. 

Figure |4] shows the Hi flaring derived in a self-consistent 
way from the LAB survey for a spheroidal halo with q — 1 in 
comparison to flaring curves expected for mass models with 
various q parameters. We discuss q — 1 first. For R <: R© 
the predicted flaring increases linearly wit h distance, in agree- 
ment w ith the analytical solution given bv lOlling & Merrifield! 
(2000, Eq. 8). Comparing the model with the data, we find a 
systematic S-shaped deviation. The model curve explains the 



average trend but fails to reproduce details of the flaring de- 
rived from the data around R ~ 17 kpc and R ~ 35 kpc. 

To improve the model we explored how the flaring might 
be affected by the flatte ning parameter q. According to 
Oiling & Merrifieldl d2000l) . we expected changes in slope. 
Oblate halos, q < 1, lead to a low flaring at R ~ 17 kpc but 
fail to reproduce the observations at R ~ 35 kpc. Prolate halos 
improve considerably for q <: 8 but cannot reproduce the data 
at R ~ 17 kpc. Each halo model also leads to a different flar- 
ing curve when fitted to the observations. Since each model 
leads, iteratively, to a different result for the rotation curve, 
there would be a need for each q parameter to also have a plot 
for the flaring, as derived from the data. For reasons of con- 
ciseness we refrain from showing all. Changes are in any case 
less significant than the changes in the model's flaring curves 
plotted in Fig.0] 
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Fig. 4. Self-consistent Hi flaring for isothermal dark matter 
spheroids. The thick red line gives the mean extension of the 
Hi layer that should be comparable to the solid black line, 
which represents the flaring curve expected for the model with 
q — 1. In addition, we include flaring curves for oblate halos 
with axial ratios c/a = 1/2 to 1/8 and also for prolate halos 
with axial ratios c/a — 2/1 to 64/1. 



Previous determinations of the flattenin g parameter q from 
H i flaring led to contradictory results, lolling & Merri field 
d2000i l200ll) consider the H 1 flaring at R = 2 R© and ad- 
vocate q = 0.8. They conclude that significantly flattened ha- 
los are possible only in the case that the dista nce of the Sun 



from t he Galactic center is smaller than 6.8 kpc. lNarayan et al 
(2005) used data in the range 8 ^ R £ 24 ^ kpc to find that 
isothermal spheroids do not fit the observed flaring data par- 
ticularly well. Using only the flaring at R — 2 R© , we could 
conclude that the halo must be rather oblate, q ~ 0.1, but Fig. [4] 
shows that the halo model needs to become incre asingly prolate 
with d istance. This reinforces the conclusion of iNarayan et al.l 
2005h that a halo described by a single constant q parameter 



does not match the observations if one considers a broad range 
of galactocentric distances. 



P. M. W. Kalberla et al.: Gravity and the Milky Way H i gas distribution 



7 



3.2.3. Flaring caused by a dark matter disk 



3.2.4. Adding a dark matter ring 



We tested alternatives to standard spheroidal halos and found 
that it is much easier to match a model to the observations if 
one considers dark matter associated with the Galactic disk. 
An example for the derived flaring curves is given in Fig. [5] 
with wha t we ex pect from our initial dark matter disk model 
(DMD) (PaperJ). For comparison we plot the mean flaring de- 
termined from the derived n(R, z, <p) distribution for the north- 
ern and southern parts of the Milky Way disk. 
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Fig. 5. Derived H i flaring for the initial dark matter disk model. 
We include the flar ing curv e as expected from the initial dark 
matter disk model (Paper I, black dashed) and after adding a 
dark matter ring to the model with parameters as given by 
de Boer et all d2005l) (black dotted). 



Comparing the flaring derived from observations with the 
model, we generally find reasonable agreement between the 
two, comparable to the highly prolate spheroidal mass model 
(q ~ 8) from Fig. |U but with a better match for 12 < R < 25 
kpc. However, we still have significant S-shaped deviations. We 
investigated whether these deviations could be explained by re- 
duction problems. Each of our attempts to explain the discrep- 
ancies with problems during the data processing failed. Even 
trying deliberately to bias the data reduction led to no signif- 
icant changes for the resulting H i distribution. We conclude 
that the deviations between model and data are highly signifi- 
cant and not caused by a systematic bias. 

The depression in the flaring curve at 12 < R < 25 
kpc has been seen befo re, see e.g. the recent papers by 
Nakanishi & Sofud d2003l) and lLevine et all d2006<J). bu t the ef - 
fect shows up as most pronounced in Fig. 23 of lVoskesI fl999). 
A low flaring can be explained either by an extraordinary low 
velocity dispersion of the Hi gas or by a local enhancement 
of the mass distribution. The turbulent velocity dispersion in 
the interstella r medium (ISM) appears to be indep endent of 



distance (Blitz & Spergell 11991b iBurton et al., 1992). We find 



None of the well-established dark matter models has an en- 
hancement of the mass distribution at a distinct radial distance 
from the center, something our models seem to require to re- 
duce or even eliminate the S-shape deviation shown in Fig. [5] 
However, a possible scenario could be a ring-like mass con cen- 
tration caused by a merger ev ent; see Hayashi et al. ( 2003b and 
Penarrubia et al.1 (120051 120061) . Several groups claimed to have 
detected stellar streams in the Milky Way, right at the region of 
interes t: |Newberg et al] d2002h;llbata et al] d2003l):lYannv et al] 
d2003b: [Martin et al.1 d2005l) : iMartmez-Delgado "et al] d2005l) : 
Bellazzini et alJ (120061) . Ot hers explain th e over-densiti es as 
due to the Galact i c war p: Momany et al. (2004, 20061) and 
Lopez-Corredoiral (120061) . The stellar mass of th e giant stel 
lar structure is estimated to be 2 10 8 M Q -10 9 M Q dlbata et al. 



2003b . It is, however, readily shown that such a mass is insuffi- 
cient for explaining the low flaring at R ~ 18 kpc. 
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Fig. 6. Observed mean HWHM flaring (thick red line), apply- 
ing the rotation curve of our best-fit model. The dotted blue 
and green lines give the flaring in the northern and southern 
hemispheres, respectively. The best-fit flaring model, contain- 
ing dark matter disk and ring, is plotted with a black dotted 
line. The flaring curve from the initial dark matter disk model 
(without ring) is shown for comparison (black dashed). 



A ring with a mass content well above the estimates for 
the stellar mass contribu tion is needed, and exactly this was re- 
cently claimed by de B oer et al. (2005) to exist. They studied 
the diffuse Galactic y-ray background observed with EGRET. 
They find excess emission close to the Galactic disk and claim 
that the excess originates from ring-like dark matter distribu- 
tions. One of the rings is located at R ~ 14 kpc. 



therefore no reason to assume that the H i gas is special at dis- 
tances 15 5 R 5 20 kpc. 



To test the Ide Boer et al.1 (12005b proposal, we include the 
ring at R ~ 14 kpc in our model (Pig. [5j. The additional dark 
matter ring indeed provides a step in the right direction, but 
obviously the predicted flaring is now too low. We varied the 
position of the ring, the radial dispersion, and its total mass. All 
calculations included a self-consistent solution of the Poisson- 
Boltzmann equation, a derivation of the Galactic velocity field, 
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and finally a redetermination of the 3-D volume density distri- 
bution n(R, z, 4>) and the resulting flaring. 

The best result is shown in Fig. [6] Note that not only the 
model flaring curve, but also the flaring derived from the data 
has changed in comparison to Figs. [3] to [5] The model repro- 
duces the observations for a ring-like mass distribution 



n(R, z) - n r e 



(R-R r )/(2a- 2 r )-\z/Zr\ 



(3) 



with n r = 0.005 cm 4 , R r = 17.5 k pc, cr r — 5 kpc, an d 
z r — 1.7 kpc. This relation was given bv lde Boer et al. ( 2005 ). 
but all parameters, except z r = 1 .7 kpc, have been redetermined 
and modified by us for an axisymmetric model. Discrepancies 



betwee n our solution and parameters given by Ide Boer et al 
(2005) may be explainable, at least partly, by uncertainties and 
systematical deviations of the ring-like mass distribution from 
axisymmetry. We determined here the mean properties of the 
ring alone. In Sects. 8 & 9 we discuss fluctuatio ns in more de- 
tail. The method used by Ide Boer et al. ( 2005k is most sensi- 
tive in direction to the anti-center, a region which is seriously 
affected by velocity crowding, therefore not accessible to our 
analysis. 

The resulting ring has a total mass of 2.3 10 10 Mq. For 
the remaining part of the dark matter m odel we kept the model 
parameters as described in Table 1 of Paper II with only two 
changes. The local dark matter disk density had to be adjusted 
from n = 0.69 cm -3 to n — 0.75 cm 3 , and we adopted a radial 
scale length Hr/2 - 3.75 kpc for a better fit of the Hi gas 
surface density. 

3.2.5. Spheroidal halo with dark matter ring 

Next we tested whether the inclusion of a dark matter ring can 
also lead spheroids to fit better. We consider three cases, q =2, 
4, and 64. Figure Q shows that a prolate spheroid with q — 4 
represents a good fit, but observations cover the range q — 64 
for the northern and q = 2 for the southern part of the Milky 
Way. The ring parameters used are identical to those discussed 
in Sect. 3.2.4, except that we need only 80% of the total ring 
mass in the case of a spheroidal halo. Solutions with differ- 
ent flattening parameters q are nearly identical for R H 25 
kpc. This emphasizes once more the need to get constraints 
at large distances. Milky-Way mass models are ambiguous if 
a limited data ba s e is used, a point stressed previously by 
Dehnen & Binneyl dl998h . 



Comparing Fig. [7] with Fig.|6]formally results in an accept- 
able solution for q = 4. Further improvements to the model 
(e.g. at R ~ 15 kpc and R ~ 22 kpc) would be possible by 
splitting the dark matter ring into several lumps. We did not 
explore them for two reasons, i) Strongly prolate spheroids 
are rather unexpected from cold dark matter simulations (e.g. 
iDubinski & Carlberg , 1991 ), and ii) the surface densities at R ~ 
Rq are unacceptable for all of our spheroidal models. In par- 
ticular, highly prolate spheroids show the most significant de- 
viations from HIPPARCOS results, see to Sect. 6 for further 
discussion. 
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Fig. 7. Derived Hi flaring for an isothermal dark matter 
spheroid with a dark matter ring. The symbols are the same as 
in Fig. [6] For comparison, flaring models are given for q — 2, 
4, and 64. 



3.2.6. NFW halo with ring 



We also tried to fit a ACDM-based NFW model (Navarr o et al. 
1996 ). using parameters provided by Klvpin et alj ( 20021) . The 
result is shown in Fig.|8]and fits the observations as well as the 
results presented in Figs. [6] & [7] and might be improved by as- 
suming a lumpy ring structure. For a favored halo concentration 
parameter C — 12, we needed, however, to modify the central 
mass density to keep v rot = 220 km s _1 , resulting in a total halo 
mass of 5 10 11 Mq. This is in clear conflict with our constraint 
that the total mass of the Milky Way within R S> 350 kpc should 
be close to 2 10 12 Mq. M oreover, the halo m ass is only 50% of 
the solution demanded by Klypin et al.l ( 2002 ), likewise in con- 
flict with the ACDM constraints and therefore not acceptable. 
We were unable to match the flaring at large distances exactly 
with the observations (Fig. [8). It would be necessary to increase 
the halo mass, but this would imply at the same time that the lo- 
cal circular speed must significantly exceed v rot = 220 kms -1 . 
For the NFW model considered, there is too much mass in the 
inner part of the Galaxy. This is the well-known cusp problem. 

4. Galactic rotation 

Each of the mass distributions discussed in the previous sec- 
tion causes a different rotation curve. This was taken into ac- 
count when we derived self-consistent flaring curves for these 
models. Here we restrict the discussion to the dark matter disk 
model. 

The 3-D rotation curves for this model are shown in 
Fig. |9l We plot the initial dark matter disk, the proposal by 



de Boer et al.l (12005b . and the final best-fit. In each case we 
show the circular streaming velocities at \z\ — to 5 kpc 
(from top to bottom). Including a dark matter ring in our model 
causes the initially falling rotation curve of the dark matter disk 
(DMD) model to rise significantly at R ~ 18 kpc. Adjusting the 
mass density of the ring to match the Hi flaring brings the 
rotation curve down. Our best-fit mass model leads to a sur- 
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line 



8. Derived H i flaring for the NFW model. The solid black 
gives the model. See Fig.|7]for further explanation. 



prisingly simple result. The rotation law at z — kpc is almost 
flat over a range 5 £ R ^ 27 kpc. 



(e.g. Blitz 


1979 


Schneider & Terzianlj 19831 Fich et al. 


1989; 


Merrifiek 


1992 


Brand & Blitzlll993l:lHonma & Sofue 


1997) 


which instead result in the high velocities of the de Boer et al. 



(2005) model. However, the discrepancy can be solved im- 
mediately if one takes distance uncertainties into account. 
The tracers used by various authors to measure the Galactic 
rotation curve for R S 11 kpc have large distance errors. 



Binn ev & DehnenI d 1997b assumed that most of these objects 
are located within a ring at R ~ 14 kpc and derived for this 
case a flat rotation curve. 
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Fig. 9. Rotation curves for the models displayed in Figs. [5] & 
[6] The black dashed lines give the rotation curves for the ini- 
tial dark matter disk, and the dotted lines the curves for the 
dark m atter disk including a ring according to Ide Boer et al. 
(2005). The thick red lines show the rotation curve for the best- 
fit model. To the rotation at z = kpc (top thick curves) we add 
in each case the circular streaming velocities at \z\ = 1 to 5 kpc 
(from top to bottom). 



4.1. The co-rotation problem 

In Fig. [10] we display the volume density distribution at a con- 
stant galactocentric distance R = 10 kpc, derived for the stan- 
dard model with a constant rotation velocity v rot = 220 km s _1 . 
FigureQT]displays the distribution at the same distance, but for 
our best-fit model. Discrepancies are obvious, because they are 
caused by the different rotation curves. The velocities at z = 
kpc are identical for both cases (see Fig. [9). Correspondingly, 
the densities at z — kpc are the same. The standard model 
(Fig. [TO! has a cylindrical rotation, but for the best-fit model 
(Fig. [TTJi v rot (R,z) lags according to Eq. 1 for increasing z dis- 
tances increasingly behind the disk. 

We emphasize that Figs.[10]&[TT|have been derived by run- 
ning identical procedures, the only difference being the 3-D ve- 
locity field. Comparing distinct features off the plane in the 3- 
D volume density distribution n(R, z, 4>) for both rotation laws, 
we find that structures, most noticeabley cloud complexes, are 
shifted to larger R and \z\ distances if one assumes a cylindrical 
rotation. An example, the HVC complex R, was discussed in 
Sect. 3.1. 

The average H i volume density of the disk < n(R, z) > falls 
off exponentially with radius R and height \z\. Shifting distinct 
features to larger R and \z\ distances therefore leads on average 
to significantly higher densities for \z\ > kpc. This effect is 
obvious if one compares the volume density distribution at a 
constant distance R as derived for cylindrical rotation with that 
derived from a lagging rotation model. The flaring is larger if 
one assumes cylindrical rotation (Figs.[3]and|6]l. Discrepancies 
are most significant close to the solar circle, but Figs. fT2l and 
[l3]derived at R = 12 kpc demonstrate that even larger distances 
may be affected. 

For our analysis we chose the lagging rotation according to 
Eq.[TJsince it considerably improves the internal consistency of 
our results in comparison to a cylindrical rotation. We neglect 
buoyancy and any depe ndency of gas pressure terms Vp/p on 
radius dKalberlal 120031 S ect. 5.5.1) . These second-order cor- 
rections, discussed also by Benjamin! ( 2002 ), may get important 
for extra-planar gas layers but are not expected to be significant 
for the disk gas. For example, at R = 20 kpc, we expect a 10% 
effect from the Vp/p term only for z ^ 14 kpc. This z distance 
is an order of magnitude off from the scale height of the disk 
gas so we ignore this correction for computational reasons. 

Our preference for a lagging rotation with increasing z dis- 
tance is supported by observations that probe extra-planar gas, 
located on top of the standard disk gas we consider here. The 
average rotational streaming velocity in the Milky Way, as de- 
rived from interstellar absorption lines at z distances up to a 
few k pc, is found to lag behin d the disk ddeBoer & Savage , 
198 31 Isavage. Sembach. &~Lulll997h . lRandl dl997l l2000h and 



Heald et al 



(2006aH)" find evidence of lagging ionized gas in 
the halos of NGC 891 & NGC 5775. Very sensitive Hi m aps 



of NGC 891 dSwaters et all J 19971 IFraternali et all 120051) and 



NGC 2403 (IFraternali etall l200ll |2004 similarly support 



lagging extra-planar H i gas phase. 
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Fig. 10. Volume density distribution of the standard model, as- Fig. 12. Volume density distribution derived for/? = 12 kpc for 
suming cylindrical rotation, derived at R = 10 kpc for azimuth the standard model with cylindrical rotation. 
180° > cf> > -180° and -15 < z < 20 kpc. Densities affected 
by local emission with |vlsrI < 10 km s _1 are flagged dark. The 
white lines indicate latitudes \b\ = 30°±0.5°. Only densities be- 
low n < 0.5 cirT 3 are displayed. We use the same logarithmic 
transfer function as in Figs.Q]&|2] 
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Fig. 13. Volume density distribution derived foxR = 12 kpc for 
the best-fit model with a lagging rotation law. 



Fig. 11. Volume density distribution derived at R = 10 kpc for 
the best-fit model with a lagging rotation according Eq. 1 . For 
more details see Fig.fTOl 



5. Consistency check - available flaring data 

5.1. Gaseous flaring inside the solar circle 

Our approach does not allow an accurate determination of the 
flaring inside the Solar circle, but we may compare our de- 
rived flaring with more accurat e data. F igure [T4lshows the flar- 
ing for the H i gas derived by ICelnik et alj dl979l) (•) and by 
Malhotral dl995l) (*). The results from these groups are found 
to be largely discrepant. 

The scale height, derived by us from the model, is the scale 
height for a mix of the CNM and WNM. Both phases have 
significantly different scale heights. ICelnik et al. dl979h have 



shown that the Hi volume density distribution perpendicu- 
lar to the Galactic plane can be described by a sech ff (z/zo) 
law. This relation can be approximated by a Gaussian distri- 
bution exp[-z 2 /(2z^)] close to t he plane and at large distances 
by an exponential, exp(-z/zA iBahca ll (1984) demonstrated 
that such an approximation also applies for a multicomponent 
model if the gravitational potential can be approximated by 
a single dominating mass component. This approximation re- 
main s valid in the case of a more general m ass distribution 
(see IPaper j. Sect. 5.2.1). ICelnik et all d 1 979b fitted a general 
sech a (z/zo) volume density law, solving for z g and z e . This is 
the most accurate and unbiased way to determine the H i scale 
height. The determination of the first moment by us, excluding 
H i gas ex tending to large z distances, is also unbiased. The ap- 
proach by iMalhotral d 1 995b is quite different. She assumed that 
the Hi layer can be fitted by a single Gaussian component, 
represented only by the WNM (Fig. [Pil l. Taking the different 
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Fig. 14. F l aring of the H i gas in the inner Ga laxy. The data, 
Malhotra d 19951) . (*), and ICelnik et all d 19791) . (•), represent 
the observed HWHM of H i clouds. The dashed black line rep- 
resents the HWHM as derived for a multiphase H i medium, 
the green dashed line is for the CNM, and the red solid line 
represents the WNM. 




approaches into account, the flaring of the H i gas within the 
Solar circle is consistent with our model. 



5.2. Molecular gas 



molecular gas flaring 
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Fig. 15. Flaring of the molecular disk. The solid line represents 
the flaring for the fina l model. The dat a poi nts give CO scale 
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In Fig. [15] we compare scale heights for the molecular gas 
phase from the best-fit dark matter disk model with observed 
scale heights. In the inner Galaxy two inde pendent analyse s 
are available from lBronfman et al.l ( 1988) and Malho"tral dl994 . 
Both sets of values appear to be slightly offset from each other. 
This may be explainable by details of the data analysis but is 
certainly unaffected by complications from a multi -phase corn- 



Galaxy we used scale heights as derived by Wo uterloot et al 



(1990) br a verification of the model. In general we find good 
agreement between model and data except for the outermost 
data points. At distances R S 15 kpc, the observations may be 
limited in sensitivity. 



6. Consistency check 
densities 



local mass surface 



Among the various observational constraints that may be used 
for determining the dark matter properties, the surface den- 
sity Si i within 1 .1 kpc is the measure wi th the lowest rela- 
tive uncertainties (Kuij ken & Gilmorelll99lh . It therefore pro- 
vides one of the most important constraints for models. Recent 
surface-density determinations based on HIPPARCOS have 

ad ded further constraint s . 

iKuijken & Gilmord dl99l[) determined Eji = (71 + 
6) M0pc~ 2 using K dwarfs. iHolmberg & FlynrJ (2004) ana- 
lyzed HIPPARCOS K giants and find E u = (74 + 6) M pc~ 2 . 
Our best-fit model provides a consis tent solution for Ey = 
(79 ± 2) M pc- 2 . Within 800 pc IHolmberg & Flvnnl d2004 



determine a surface density of Eg = (65 ± 6) M0pc~ 2 , we 
get £,x = (66.9 ± 2) M <T>pc~ 2 . For distances up to 350 pc, 
Holmberg & Flvnnl d2004 determine £35 = 41 M Q pc~ 2 , and 
Korchagin et all d2003l) find £.35 = (42 + 6) M pc~ 2 . Our result 
is S.35 = (42.7 + 2) M pc~ 2 . 

For z — pc, we find an Port limit of pp) _ . 097 M mpc 3 



in excellent agreement with IHolmberg & Flym] d2000h . who 
obtained the same value, and only marginally dev i ating from 
p© = 0.10 M©pc~ 3 derived bv lKuiiken & Gilmorel |l989) and 
PO = (0.105 ± 0.005) M pc- 3 bv lKorchadn et all d2003l) The 
total surface density for all disk material (stars and gas) in 
our best-fit model is 52.5 Mr^pc 2 , almost identical with th e 
model estimate 52.8 M©pc~ 2 by Holm berg & Flvnnl d2004 . 
Summarizing, the local volume and surf ace densities are in 
excell ent agreement with constraints by IKuijken & Gilmore 
dl99ll) and with all constraints derived from HIPPARCOS. 
The uncertainties given for our model are based on typical 
parameter changes corresponding to modifications of the ring 
parameters within an acceptable range. Our initial dark mat- 
ter disk model without ring had Eu = 72.5 Mqpc~ 2 , but 
addi ng the dark matter ring parameters as initially proposed 

, u = 88.3 M pc~ 2 , too 



by Ide Boer et al.l d2005l) resulted in E 



high in comparison to the HIPPARCOS result of E u = (74 + 
6) M0pc~ 2 . This model got acceptable only after tuning with 
respect to the observed Hi flaring (Sect. 3.2.4). 

All of the conventional halo-mass models, discussed in 
Sects. 3.2.2, 3.2.5 and 3.2.6 have surface densities E; 1 that de- 
viate by at least 30 % from Eu _ (74 ± 6) M p>pc~ 2 , derived 
from HIPPARCOS dHolmberg & Flvnnl |2004 . Extreme pro- 
late spheroids (Figs.|4]and|7]l show the worst deviations. 

7. The Galactic warp 

It is well known since Westerhoutl ( 19571) that the Milky Way 



H 1 disk is warped. We need to take this into account. The grav- 
itational potential for our model is calculated for a flat system, 
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but we apply a geometrical transformation afterwards to warp 
the system. We describe W(R, 0) as the deviation of the center 
of the gas layer from z — by 

W(R,cf>) = W (R)+Wi(R)sm(<p-<p 1 (R))+W 2 (R)sm(2<p-(f>2(R)),(^) 
where W, and 0, the amplitudes and phases of mode i. 
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Fig. 16. Amplitudes for the warp modes to 2 derived for the 
best-fit model. Mode 2 is insignificant for R < 15.5 kpc. 
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Fig. 17. Phases for the warp modes 1 and 2 derived for the best- 
fit model. Mode 2 is insignificant for R < 15.5 kpc. 

To derive the warp we averaged the observed densities 
within rings having a width of AR = 1 kpc for R = 10 to 40 kpc. 
For each annulus we fit the modes to 2 simultaneously using 
a least squ are sine wave fit b ased on the Marquard/Levenberg 
algorithm ( Press et allll986t) . Mode 2 is found to be insignifi- 
cant for R < 16 kpc. The fit in this range is therefore without 
mode 2. We found that azimuth angles 80° $ 4> $ 130° for 
R 'k, 16 kpc show significant perturbations and ex clude this 



range (region X as discussed bv lLevine et al .. 2006a). 



The formal errors for our fitted parameters are quite low, a 
few percent for the amplitudes, and typically less than a degree 



for the phases. Only at large distances do the errors increase 
for mode 2. At R ~ 40 kpc we find a 10% amplitude error, the 
typical mode 2 phase error is 3°. 

Figure[l6]displays the ampli tudes for the basic w arp modes. 
We find a good agreement with lLevine et al.1 d2006al) . In Fig. [171 
we display the warp ph ases. We followed the convention used 
bv lLevine et al.l (12006a) and plot the phases for the first maxi- 
mum close to azimuth 90°. Here we find some differences for 
the phase of mode 2, most probably due to the different rotation 
curve used by us. In Fig.[l7]we see a well-defined maximum in 
the phase of mode 1 between 15 % R % 25 kpc. This appears 
to be associated with an almost constant slope in the phase of 
mode 2. Formal errors for the phases in this range are very low, 
less than one degree. 

The warp modes are well-defined. The amplitude of mode 
rises continuously up to R ~ 40 kpc. Modes 1 and 2 both 
show a comparable trend in their amplitudes, peaking close to 
33 kpc. The phases are also quite regular. 

8. Deviations from axisymmetry 

We modeled the flaring of the Milky Way Hi gas assuming 
that the distribution of mass and gas in the Milky Way is in 
equilibrium. An axisymmetric model was used to fit the mean 
flaring. Figures [3] to [8] indicate systematic differences for az- 
imuth angles <p ^ 180° (northern sky) and <p ^ 180° (southern 
sky), and all flaring curves show in general the same trends: 
The axisymmetric model does not reproduce the flaring exactly 
but represents average properties of the mass distribution, im- 
plying a massive exponential disk with an additional massive 
ring at R ~ 17.5 kpc. Here, we intend to discuss systematical 
perturbations from axisymmetry. 

8.1. A lopsided disk 

We calculated the average flaring within 20° wide sectors. 
Figure [18] shows for the southern part of the Galactic disk that 
the sector averages are well-defined, deviating only a few per- 
cent from the mean. In the northern part (Fig. [T9l the situa- 
tion is very different. For radial distances R ^ 25 kpc and az- 
imuth 50° S> <p ^ 130°, we find a significant increase in the 
flaring. The strongest deviation from the averag e flaring is at 
90 ° ^ <j> % 110°. It has already been noted by iLevine et al 



(2006a) that this region is special. Inspecting the survey data, 
we find a highly disrupted H i distribution with some evidence 
that the H i in this range may be significantly perturbed (see 
Figs.Q]&[2|i. We exclude this sector from the average over the 
rest of the disk. 

Systematic deviations in the flaring between northern and 
southern part of the disk are explainable if either the mass dis- 
tribution is non-axisymmetric or the state of the ISM is deviant 
for a part of the disk. Explaining the asymmetries for R ^ 25 
kpc as merely due to the mass distribution leads to a dark mat- 
ter model that in the southern part has three times the mass of 
the northern hemisphere (the fits S and Nl in Figs. [T8l&|19ll. 
Correspondingly, the rotation velocities at R ^ 25 kpc should 
vary by 50 km s _1 . This is not observed so we exclude this pos- 
sibility. 
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Fig. 18. Flaring of the southern part of the Galactic disk. The 
average flaring for 205° < (f> < 350° (dashed green) is com- 
pared with the model fit (thick red). Flaring curves, derived for 
individual 20° wide sectors, show good internal consistency of 
the flaring for the southern part of the disk. 
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Fig. 19. Flaring of the northern part of the Galactic disk. The 
average flaring, excluding 50° < <p < 130°, is indicated by the 
dashed blue line, the corresponding fit (Nl) by the black line. 
Flaring curves derived for 20° wide sectors show systematical 
deviations for R > 25 kpc. The thick red line shows the fit 
(N2) to the low flaring part of the H i , for 30° < <p < 50° and 
130° <<h< 170°. 



Taking the opposite point of view that the flaring is affected 
by variations in the state of the ISM, we assume for the moment 
that the flaring for 50° ^ ^ 130° is "atypical" and may be 
disregarded (a more detailed discussion is given in Sect. 9). For 
the low flaring part of the northern disk we find that the dark 
matter disk has to be somewhat more extended (radial scale 
length Hr/2 = 8.75 kpc) than the disk fitted to the mean flaring 
used in the previous sections (Hr/2 = 7.5 kpc). For the south- 
ern part of the disk (205° < <p < 350°), we find almost identical 
parameters, Hr/2 = 8.5 kpc. Fitting the low flaring sectors only 
(model N2 & S) results in a dark matter disk that is, within the 
errors, axisymmetric. The total mass of the dark matter disk 



is affected by the choice of the model. For best-fit model as 
displayed in Fig. [6] we obtain a mass Momd = 1-8 10 11 Mo- 
Fitting only the low flaring part of the H i disk (fits N2 & S in 
Figs.0J&[T9]l leads to M DMD = 2.4 10 11 M Q . 

Major north/south differences between the models remain 
and are easily visible in Figs. [18] & [19] at R ~ 25 kpc. This is 
caused by variations in the dark matter ring radius with best- 
fit results of R - 13 kpc in the north and R = 18.5 kpc in 
the south. Similar variations between 15 and 20 k pc for the 
observed stellar ring structure have been noted by Ibata et al. 
(2003). The total mass of the dark matter ring is M — 2.2 to 
2.8 10 10 M . Fits Nl and S imply further that the southern 
part of the dark matter ring is three times as massive as the 
northern part. Our finding is consistent with know n stellar over- 



densities in the south (e.g. Bell azzini et al. . 20061) . suggesting 



an association between dark matter ring and stellar distribu- 
tion. The ring, however, cannot be associated with the Canis 
Major over-density. The radius R = 18.5 kpc of the dark matter 
ring in this direction implies that this object is located behind 
Canis Major and ins tead associated with "the One Ring" dis - 
cussed previously by Newber g et alJ (200i) :llbataetaljl2003h : 
Rocha-Pinto et all d2003l) and lYannv et alJ d2003b . 

North/south differences in the mass distribution, consider- 
ing the mass distribution from dark matter ring and disk allto- 
gether (fits N2 & S), lead to a lopsided dark matter distribution 
that is most significant at R ~ 25 kpc. Our fits imply that the 
rotational velocities at this distance deviate by ~ 15 kms~', in 
good agreement with the empirical firs t -order e picyclic stream- 
line correction derived by Levine et al. ( 2006ah . We applied this 
correction as described in Sect. 2.1. 



8.2. Spiral structures 

For a better understanding of the nature of the large-scale 
fluctuations, we display in Fig. [20] the relative deviations of 
the observed flaring from the axisymmetric model defined as 
HWHMofcj/HWHMmode/. The north-south asymmetry is obvi- 
ous, but we also find ring-like perturbations and spiral arm fea- 
tures. Comparin g Fig. l20l with perturb ations in gas thickness 



as determined bv lLevine et a l. (2006b, Fig. 2) shows excellent 



agreement, although very different methods have been used. 
This applies to the filtering of the data base for the separation 
of disk emission, to the determination of the scale heights, and 
in particular to the derivation of perturbations in scale height. 
Levin e et al. (2006b) use the un-sharp masking method, while 
we compared the data with a model. The only agreement be- 
tween both methods is that similar rotation curves were applied 
to conver ting the observed spec tra to a 3-D volume density dis- 
tribution. iLevine et al. I d2006allbh assume a flat rotation curve, 
v = 220 km s _1 , while we derived the galactic rotation (Fig. [9]l 
by fitting a mass model, also ending up with a flat rotation law 
for z = kpc. 

First of all we need to discuss the north-south asymmetry 
of the flaring. This begins at R ^ 15 kpc, about the same as the 
distance at which the warp causes a significant bending. The 
four-arm spiral pattern according to ILevine et al.l (12006bl) . in- 
dicated in Fig. [20] also shows a marked asymmetry with arms 
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wide open to the north. Structures like these suggest tidal inter- 
action. The Magellanic system might have caused a dark mat- 
ter wake, enhanc ing the mass distribution towards the south 
(Weinber g[ ll998h . The tidal interaction scenario is consistent 
with the onset of the warp and higher- mode vertical pertur- 
bation s at larger distanc e s (|Saha&Jog , 20061: Levine et al. . 
2006 at IWeinberg & BlitzL 120061) . Large-scale deviations from 
axisymmetry, visible in Fig. [20] are fully consistent with such 
a scenario. 
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Fig. 20. Deviations HWHM ofo /HWHM morfe/ for the observed 
flaring relativ e to the best-fit mode l. The location of spiral arms 
according to iLevine et al. (2006b) is indicated. The circles in- 
dicate galactocentric radii of 10, 20, 30, and 40 kpc. 



We used the s piral-arm parameters as determined by 
Levine et al. (2006b, Table 1) to verify whether the depressions 
in scale height along the spiral arms may be caused by gravi- 
tational forces from the stellar population. For the arm region 
we set a constant width of A(f> = 15° in azimuth, as indicated 
in Fig. [20] This is compared with an inter-arm region with a 
width of 17° 5 on each side of the arm. Perturbations in direc- 
tion to the Galactic center or anti-center, caused by spurious 
emission, are excluded as indicated. Figure [2T| plots the ra- 
tio HWHM,„ /e ,_ nrm /HWHM nrm determined as an average over 
all four arm s and a lso for arm 3 alone. Arm 3 according to 
Levine et al. d2006bl Fig. 3) was chosen since it is particular 
well-defined. It runs from R ~ 12 kpc, 4> ~ 180° to R ~ 40 kpc, 



4> ~ 340°. 

We used our model to calculate the flaring perturbations 
relative to an axisymmetric model, caused by enhancements of 
the stellar volume density within the arms by a factor of 3 or 6 
over the axisymmetric model. A factor of 3 would be expe cted 
for well-defined spiral arms, as in M51 dRix & Riekelll993l) . A 
factor of 6 may be considered as an upper limit in the case of a 
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Fig. 21. Flaring contrast HWHM,„ rer _ arm /HWHM arm between 
inter-arm and arm regions derived for the arms indicated in Fig. 
[20l The black lines represent model calculations assuming that 
the flaring is affected by gravitational forces from stellar spiral 
arms. 



four-arm spiral if the stars concentrate in arms with A<p = 15° 
in azimuth. 

Comparing the observed contrast with the model, we find 
that the observed amplitudes of the stellar spiral arms in the 
outer parts of the Milky Way are stronger than expected. This 
finding, however, is not solid enough to conclude that the spi- 
ral structure in the H i flaring must be caused by an additional 
mass component. Gas, according to our model, contributes only 
13% to the baryonic surface density of the outer disk. The 
self gravity of the gas is therefore insignificant. The model pa- 
rameters used for the stellar disk, in particular for the expo- 
nential radial scale le ngth of the thick disk (Hr = 7.5 kpc, 
Chiba & Beersl 120001) . are highly uncertain. The surface den- 
sities for the outer stellar disk may be underestimated in our 
model. We conclude that spiral arms influence the mass distri- 
bution out to R ~ 30 kpc, possibly even out to R ~ 35 kpc if 
we consider only arm 3. For comparison. ILevine et al . (2006b) 
find evidence of a spiral structure out to R ~ 25 kpc. The pres- 
ence of stellar spiral arms appears to be noticeable out to large 
distances, approximately 2/?25- 

Figure [20] shows some ring-like fluctuations and Fig. [2T] 
gives evidence of radial variations of the flaring contrast. A 
close inspection of Figs. [6] [18] and [19] also shows some os- 
cillations if one compares the model flaring with the observa- 
tions. Obviously the mean radius of the ring (R ~ 17.5 kpc), 
as used in our model, is only a crude approximation of a more 
complicated pattern of radial substructure s in the Milky- Way 



mass, as expected for merger events (e.g. iMartin et all 12005 



Fig. 11). The tidal interaction scenario is consistent with the ob- 
served radial variations of spiral arm amplitudes that may arise 
from interference of a pre-existing spira l pattern with tidally- 
induced spiral arms dRix & RiekeL Il993l) . Radial variations in 
the arm amplitudes in our Galaxy may be comparable to varia- 
tions found in M5 1 . An alternative interpretation of radial fluc- 
tuations are ring-like mass concentrations caused by mergers. 
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Figure [2J] shows that these also influence the arm/inter-arm 
contrast. 

In summary, we conclude that our model is able to repre- 
sent the global properties of the Milky-Way mass distribution. 
Even substructures like spiral arms and radial fluctuations are 
explicable. 

9. Uncertainties 

In the previous sections we have claimed that a dark matter 
component associated with the Galactic disk can fit the ob- 
served H i flaring best, but at least we have found no easy way 
to match a standard model to the flaring data. The best conven- 
tional model, a highly prolate spheroid (2 $ q iS 64), appears 
unlikely since this is not what we expect from ACDM models. 
Could it be that the analysis is systematically biased? 

As discussed in the previous section, the gravitational field 
of the Milky Way must have significant non-axisymmetric con- 
tributions. Accordingly, the Galactic Hi disk cannot be in an 
axisymmetric equilibrium state. In addition, stellar activities 
may lead to a fountain flow and mass may be accreted from 
outside ( van Woerden et all 120041) . As part of our analysis we 
estimate the amount of gas in non-circular motion for each po- 
sition in the Galactic plane, thus either representing out- or in- 
flow. We exclude this gas from our analysis, on average 3% to 
5% of the total H i. The remaining uncertainties are most proba- 
bly below a percent level, too low to bias derived flaring curves 
significantly. 

Most critical for deriving halo properties is the modeling of 
the Hi scale height. The conventional approach is to assume 
a singl e-phase medium with a velocity dispersion of a ~ 9.2 



km s" 1 tolling & Merrifieldl.l2000l.l200 it Naravan et aill2005l) 
We considered a more sophisticated two-phase me dium with 
two d ifferent characteristic scale heights dDickey & Lockmanl 
199(1 defined by the turbulent motions within the Hi gas . 
Using the Leiden/Dwingeloo Survey dHartmann & Burtoni 



ing 

1997), Kalberla & Kerp ( 1998) determined the average expo- 



winge 
srp (1? 



nential scale height of the CNM to 150 pc and the scale height 
of the WNM to 400 pc (see also Table 1 in lPaperJ) . The total 
scale height of the H i gas layer is in general given by the mix 
of both phases and therefore depends on the properties and the 
composition of the two phases throughout the disk. 

The conditions for the existe nce of a two-phase medium 
have recently been discussed by IWolfire et al.l (I2003h . These 
authors conclude that the CNM and WNM phases must coexist 
over most of the Galactic disk, at least out to radii of 16-18 
kpc. Basically, such a two-phase medium can exist as long as 
the turbulent pressure of the WNM is high enough to drive parts 
of the H i gas into a co l d den se phase. We repeated the calcu- 
lations by IWolfire et al.l d2003l) of the parameters of the H i gas 
distribution. We find, for large radial distances, a 5 0% higher 
turbulent pressure than used bv lWolfire et al.l d2003l) . This cer- 
tainly supports a two-phase medium at large distanc es. From a 
Gaussian analysis of HVCs. IKalberla & Hau d (2006) find clear 
evidence of a two-phase medium in complex H (R ~ 30 kpc), 
the Magellanic stream, and the leading arm (R ~ 50 kpc). In ad- 
dition they argue for a high pressure to the hot halo gas phase. 
We conclude that there is strong evidence that the H i gas has 



a multiphase composition throughout most of the Milky Way 
disk. Any reservations that the C NM must be absent beyond 
the e xtent of the stellar disk (e.g. lolling & Merrifieldl [2000, 
2001) appear to be outdated. 

The extended strongly flared region, visible in Fig. [20] at 
R k 25 kpc and 50° ^ <p S= 130°, deviates significantly from 
the average Hi disk gas. This region has particularly low sur- 
face and volume densities and it is strongly warped. It is quite 
possible that the H i gas is out of equilibrium there and that the 
CNM phase may be under-abundant. Detailed investigations 
would be necessary to verify this, but they are clearly beyond 
the scope of this investigation. 

Figure [14] displays scale heights of our model individu- 
ally for the CNM and WNM phases, but also for the two- 
phase mix. Note that the scale height of the two-phase medium 
is only slightly larger than the scale height of the CNM. A 
WNM gas phase without any cold Hi gas would result in 
a scale height that is twice as big as the scale height for the 
CNM/WNM mixture. This factor of two defines an upper limit 
for an excess in flaring caused by the transition of a multi- 
phase gas into a WNM-dominated phase. Except for the range 
90° % <p % 110°, such a factor applies to the region R <: 25 
kpc and 50° S5 <p ~ 130°, described in more detail in Sect. 8.1. 
For the rest of the disk we have no indications that there are 
significant changes in the H i composition for R ^ Rq so use a 
constant column density ratio for the CNM and WNM phases 
throughout the Milky Way. 

Calcu lating the scale heights from the mass model (Eq. 7 
in Paper I) . we are faced with the problem that individual gas 
components are not only supported by their own turbulence, 
but are also affected by the p r essure of ma gnetic fields and cos- 
mic rays. Following Parker ( 1966 , 19691) . the support can be 
described as an increase in the velocity dispersion by a factor 
yl + a + ji, where a stands for the increase due to the magnetic 
field and j3 due to the cosmic-ray component. A simultaneous 
determination of the 3-D distributions of gas, magnetic fields, 
and cosmic rays is needed to derive a and /?. F rom observations , 
the distribution of H i, th e diffuse ionized g as ( Reynolds! 1997 ). 
and the 10 615 K pl asma dPietz et all 19981) . as well as the syn- 
chrotron radiation ( Haslam et all 1982b and the y-ray emission 



(Fic htel et all 1 19941) . are ava ilable for this purpose. B ased on a 



correlation of these datasets. [Kalberla & Kerpl l 1998 ) conclude 
that, for the CNM and the WNM gas phases, only a limited 
support is expected from a turbulent magnetic field component 
close to the disk with a ~ 1/3. Accordingly, this leads to an ef- 
fective increase in the velocity dispersion by a fac tor of ~ 1.15. 
For comparison. lOlling & Merrifieldld2000ll200lh . arguing that 
the pressure support can be derived from the scale heights of 
gas, estimate a value of ~ 1 .06 for the increase in the effective 
velocity dispersion due to magnetic fields and cosmic rays. 

Our algorithm allows using separate scale heights for 
the CNM and WNM, but we find no reason to vary the 
CNM/WNM ratio or the pressure support from magnetic fields 
and cosmic rays with galactocentric radius. Our solution is 
close to a single-component model with an effective constant 
velocity dispersion of 8.3 kms , somew hat low in compari- 
son to the pr eferred value of 9.2 km s -1 of Idling & Merrifield 
(2000, 2001). These authors explored the parameter space, af- 



16 



P. M. W. Kalberla et al.: Gravity and the Milky Way H i gas distribution 



fecting H i flaring for various dispersions at a single constant 
galactocentric radius R = 2 Rq. In contrast, we use a con- 
stant dispersion for a broad range of distances, and the ques- 
tion arises whether this is justified. It is very difficult to esti- 
mate how far the isothermal assumption may lead to possible 
biases. However, we may turn the question around and ask how 
far we would need to modify our ISM model assumptions if we 
want to save the paradigm of the Milky Way surrounded by a 
spheroidal halo. 

For a spheroidal mass model without a ring and q — 1 (Fig. 
2]l we would need to increase the estimate for the flaring at 
R ~ 17 kpc by 40% and, at the same time, to decrease the es- 
timate at R ~ 35 kpc by 40%. Pressure support from magnetic 
fields might be missing at large distances, which would bring 
the scale down by 15%. In turn, we need to assume that veloc- 
ity dispersion of the Hi gas is overestimated by 25%. Given 
the uncertainties, such a solution cannot be excluded. The ar- 
guments are the opposite at R ~ 17 kpc. Here we need to boost 
the turbulent gas pressure by 40%. However, the question im- 
mediately arises as to what kinds of processes would be avail- 
able to feed this additional turbulence at distances 12 ^ R £ 22 
kpc. This region is largely outside the stellar disk, hence we 
expect no support from stellar activities. Figure H] shows that 
highly prolate spheroids fit best at large distances. However, 
the problem remains to explain a turbulent pressure support at 
R ^ 25 kpc. Alternatively, highly oblate models fit the observa- 
tions well up to R £ 25 kpc, but it is then a serious problem to 
explain the observed flaring in the outskirts of the Milky Way. 

Other mass model s without dark matter in the Galactic disk 
might be considered. iNaravan et al. ( 2005;|) advocate a model 
according to Eq. [2] with a power index p — 2. Such a model 
fits the Hi flaring well. The total halo mass is 2.8 10 11 Mq, 
however, much too low to be acceptable. This model has most 
of its mass within R ^ 50 kpc. Adding an extended halo to 



the lNarayan et al. (2005) model would result in a bimodal dark 



matter model, very similar to our case. 

10. Summary 

The aim of our contribution is to derive an H i volume density 
distribution n(R, Z, <p) for the Milky Way. The main driver be- 
hind this is the availability of th e LAB survey, a new all-sky 



Hi survey (IKalberla et al.U 2005). In particular, we expect that 



the low internal errors of this survey may lead to higher fidelity 
in the derived distribution. 

Previously, this survey was used bv lLevine et al. I d2006allbh . 
who assume a constant v = 220 kms -1 rotation law. In per- 
forming the Tb(1, b, v) to n(R, z, <p) conversion we notice, how- 
ever, that uncertainties in the Milky Way rotation curve are 
large enough to make the results ambiguous. We have tried to 
solve these problems by matching mass models of the Milky 
Way to the observations. The strategy is to obtain a self- 
consistent solution. We generated mass models that satisfy both 
the Poisson and the Boltzmann equation at the same time. To 
apply a model we require that predictions from the mass model 
need to agree with the derived volume density distribution. In 
particular, we used the global flaring of the H i gas layer to de- 
cide on the consistency of a model. The link between derived 



volume density distribution and flaring according to a mass 
model is the rotation curve, requiring a self-consistent mass 
distribution. Thus, our procedure is based on a closed loop. 

10.1. Properties of the mass model 

The Galaxy is embedded in an isothermal dark matter halo with 
a core radius of 35 kpc and a total mass of M = 1.8 10 12 
M0 within 350 kpc. To fit the global flaring of the H i disk, 
we need dark matter associated with the disk in the form of a 
self-gravitating exponential disk with a radial scale length of 
7.5 - 8.75 kpc and a mass M DMD = 1.8 to 2.4 10 11 Mq. To 
explain H i flaring at 1 3 £ R £ 20 kpc, we need an additional 
component. A dark matter ring at a radius 13 & R & IS. 5 kpc 
contributes M = 2.2 to 2.8 10 10 Mq. About 2/3 of the ring mass 
is located in the south at R ~ 18 kpc, while the remaining part 
of the ring mass is in the north at R ~ 13kpc. 

According to our best-fit model (the average flaring as dis- 
played in Fig. |6j, the total mass of the Milky Way disk within 
50 kpc is M = (2.9 ± .1) 10 11 Mq. However, the state of the 
ISM at R > 22 kpc is rather uncertain (Sects. 8 & 9). The 
northern part of the outer disk may be dominated by the WNM 
phase. Taking this into account by fitting the low flaring part 
of the Hi disk only (fits N2 & S in Figs. US & [19) leads to 
M ~ 3.5 10 11 M©, probably an upper limit. 

10.2. Global properties of the H i disk 

Our best-fit model results in a rotation curve that is essentially 
flat for R ^ 27 kpc and falling off gradually for larger dis- 
tances. Our derived Hi volume density distribu tion resembles 
previous resu lts obtained for a constant rotation dVoskes[ ["l999: 



Levine et al. , 2006al) . but we derived cleaner pictures of the 



Milky Way H i gas distribution and extended the investigations 
to larger distances. A significant improvement is obtained by 
dropping the assumption that the gas is barotropic, hence in 
cylindrical rotation. 

We derived a well-defined exponential H i disk with n(r) ~ 
exp(-R/3). The Hi surface density distribution is also expo- 
nential and has a radial scale length of 3.75 kpc. The Hi disk 
flares strongly. For R ~ 4 kpc the HWHM scale height of the 
H i layer is 0.06 kpc, increasing to 2.7 kpc at ~ 40 kpc. 

The H i disk shows systematic north-south asymmetries for 
R <; 15 kpc, apparently correlated with the warp. The warp has 
well-defined Fourier modes to 2 with slowly variable phases. 
The mean flaring in the northern part is up to a factor of 2 larger 
than the flaring in the south. 

11. Discussion 

The dark matter disk, propo sed by us, appears to conflict with 
results from HIPPARC OS dHolmberg & Flvnni l2000l I2004 ; 



Korchagi n et al.U2003b . The HIPPARCOS data are believed to 



imply that there is no evidence for significant amounts of dark 
matter in the disk. However, the flaw in this frequently used 
statement is that a definition of the disk is missing. Constraints 
from the stellar surface densities define only a local lower limit 
of 2 - 3 kpc for the scale height of a dark matter component 
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associated with the disk. The dark matter disk in our model has 
a local scale height of h z - 4 kpc and matches all HIPPARCOS 
constraints. Models without any dark matter within the disk are 
not consistent with the observed gaseous flaring, in particular 
for distances R ^ 22 kpc. 

A dark matter disk, or alternatively a highly prolate 
spheroid, explains the global flaring of the Galactic H i disk. 
But there remains a marked depression in the derived flaring at 
distances 15 $ R ^ 20 kpc. To explain this it is in any case 
necessary to add a massive ring. 

11.1. Lord of the rings 

A longstanding discrepancy, discussed in Sect. 1, is that the 
observed rotational velocities f or R <; conflic t with mass 
models. To solve this problem. iBinney & DehneiJ ( 1997h pro- 
pose that most of the tracers that appear to be at R ^ 11 
kpc might actually be concentrated within a ring at R ~ 14 
kpc. Stellar over-densities suggesting a ring-like concentration 
were found a few years later by several groups independently. 
Pe rhaps the most co ncise description of this feature was given 
bv llbata et al. ( 2003 ): there is a giant stellar ring at 15 55 R ^ 20 
kpc with a radial thickness of 2 kpc and a vertical scale height 
of ~ 0.75 kpc. The total stellar mass is estimated to be between 
2 10 8 M and 10 9 M . 

Such a mass, however, is insufficient for explaining the ob- 
served flaring. In dependent evidence of a massive ring was 



given recently by Ide Boer et al. (2005). EGRET observations 



show excess in the diffuse Galactic y-ray background. A possi- 
ble explanation is that the giant stellar ring is associated with a 
dark matter ring. We included this mass component in our dark 
matter disk model but noticed that the m ass of such a ring, ini- 
tially proposed by Ide Boer et alj d2005l) . needs to be reduced, 
and the position and extension also needs to be slightly modi- 
fied to explain the H i flaring. A fully satisfying self-consistent 
solution of our mass model and the average derived density 
distribution is obtained for a dark matter ring at R ~ 17.5 kpc 
with a radial extension of 5 kpc and a mass M = 2.3 10 10 
Mq (see Fig. |6j. The associated rotation curve is almost flat 
out to R 5 27 kpc. This extended dark matter ring, associ- 
ated with a somewhat less extended stellar ring and some en- 
hanc ements in H i is compl etely consistent with the proposal 
by Binnev & Dehnen ( 1997b that distance errors may mimic a 
rising rotation curve for R ^ 1 1 kpc. 



The lde Boer et alj ([2005) paper explains the EGRET excess 
emission as due to dark matter annihilation (DMA). From the 
energy spectrum of the excess radiation, they deduce a WIMP 
mass between 50 and 100 GeV. Th e significance of the EG RET 
excess emission was confirmed by Bergstr om et alj J2006), but 
these authors question that the excess can be cau sed by DMA 



They model the decay processes proposed by Ide Boer et al 



(2005) and derive that such a DMA would be associated by a 
significant antiproton flux. The expected flux is well above ob- 
servational limitatio ns, so they conclude th at the signal cannot 
be caused by DMA. Bergstr om et al. (2006) propose to explain 
the EGRET excess by enhancements in the baryonic matter dis- 
tribution within the area covered by the ring. 



As a crude order-of-magnitude estimate, the column den- 
sity of the ring along the line of sight is comparable to the to- 
tal column density of the observable H i disk in the same di- 
rection. For such baryonic matter contained within the usual 
gas phases, a huge EGRET signal would be expected. This 
is clearly incompatible with the background signal observed 
by Ide Boer et alj J2005I) . where any possible baryonic matter 
within the ring must be largely hidden from the observer by 
self-shielding effects. 

Baryonic dark matter in the form of cold, dense molecular 



gas clumps (clump uscules), associated w ith th e H i disk, was 
proposed in iPfenniger & Combes (1994) and IPfenniger et alj 



d!994l) . Such clumps are hardly detectable optically or with ra 
dio telescopes. They are self-shielding against y-ray emission 
if their diameters are suffi ciently small, e.g. ^ 10 AU for a 
clump mass of 10 3 M© dKalberla et all 1 19991: lOhishi et al 



2004). Excess y-ra y emission, correlat ed with local gas and 
Ir 



dust was found by brenier et al. (2005). These investigations 
suggest that some "dark gas" is associated with Hi gas, dis- 
tributed like an envelope around CO-traced molecular gas. The 
total mass of this "dark gas" is comparable to the mass of the 
molecular gas, negligible in comparison to the dark matter ring. 

1 1.2. Origin of the dark matter ring and disk 

Our mass model is axisymmetric. However, the derived flaring 
shows significant asymmetries (Figs.[T8ltol20l. Spiral structure 
exists out to R ~ 30 kpc, and we find strong radial variations 
and a pronounced north-south asymmetry. Typical scales for 
systematic fluctuations appear to be in the range 3 to 15 kpc. 

The Milky Way is not a unique galaxy. Rotation curves of 
other galaxies show simila r wiggles; perhaps the most impres- 
sive collection is given by Sofu e et all {l999). It appears hard 
to explain all these fluctuations as caused by an extended dark 
matter halo. The current understanding is that the halos rep- 
resent a rather smooth distribution. Despite density enhance- 
ments at the center (cusps), no substructures on kpc scale are 
expected. Clumpy ness appears in stead to be associated with 
dissipative matter ( Combes , 2002h . 

The dark matter ring may provide a key to understand- 
ing these fluctuations. The ring, located at 13 :S> R & 18.5 
kpc, has a FWHM extension of 12 kpc. Associated stellar 
streams appear to be more sharply bounded in distance, in 
FWHM typically 2 kpc to 4 k pc dRocha-Pinto et all 12003 



hWHM typically z kpc to 4 k pc qjK ocha-Finto et all IZUUJI 
iMartmez-Delgado et all 120051) . iNewberg et alJ d2002l) deter- 



mine an upper limit of 6 kpc. Simulations of the accretion of a 
dwarf onto the Galaxy show th at the most promine nt tidal arms 
are also spatially well-defined dMartin et al. , 2005 ). 

If the giant stellar ring is caused by the disruption of a 
dwarf galaxy, it appears reasonable to also associate the dark 
matter ring with the dwarf. The mass-over-light ratio for the 
progenitor galaxy would be in the range 24 to 120, with the 
dark matter distribution more extended by a factor of 4. This 



case re sembles the Sagittarius dwarf, as modeled bv llbata et al 



dl997|) and llbata & Lew is ( 1998). These authors needed a sim- 
ilar dark matter halo to stabilize the dwarf against the Galactic 
tides. 
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Helmi& White! d 19991) argue that the stellar halo may have 
been built up by disrupted satellites. Their calculations sug- 
gest that the stellar halo could consist of 300 to 500 streams. 
Assuming that these streams are associated with dark matter 
would imply a more or less smooth spatial distribution of dark 
matter debris within the Galaxy. The hypothesis that satellites 
may have conside rable amounts of dark matter is supported 
by iGilmore et al who argue that Milky Way dwarf 

spheroidals are dark matter dominated by halos extending pos- 
sibly to k pc scales with typica l mass-to-light ratios of a few 
hundreds. lHavashi et al.l d2003l) essentially come to the same 
conclusion when simulating the structural evolution of sub- 
structures. From dissipation-less simu lations on structur e for- 
mation, a roundish halo is expected dAllgood et al.Ll2006l) . The 
dark matter disk advocated by us implies a strongly anisotropic 
distribution for the progenitors of the stellar streams. The orbits 
for the majority of the accreted dwarfs must have been highly 
correlated with the plane of the Milky Way disk. 



Will man et al.l (2004) discuss the completeness of the cen- 



sus of Milky Way companions with respect to Galactic lati- 
tudes. For a random distribution, they estimate an incomplete- 
ness as large as 33%. Up to 10 satellites at low latitudes, pos- 
sibly causing 300 to 500 streams, if captured, could be waiting 
for detection or are already captured by the G alaxy. The recent 
detec tion of five new Milky Way companions (Belokurov et al., 
120071) implies that this is only a lower limit. lAbadi et al.l d2003) 
demonstrate for disk galaxies that the bulk (~ 60%) of the thick 
disk consists of the tidal debris from satellites whose orbital 
plane was coincident with the disk and whose orbits were cir- 
cularized by dynamical friction prior to full disruption. If these 
satellites contained also dark matter, one would expect that this 
material would have built up a thick disk-like dark matter com- 
ponent. The shape of the resulting dark matter disk would de- 
pend on the accretion history and on the nature of the dark mat- 
ter. 
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